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1 Introduction 

Equilibrium phase transitions, e.g. the transition at 0°C from water to an ice- 
crystal, are common phenomena in nature. Such phase transitions between a 
disordered high temperature phase and an ordered low temperature one are 
rather well understood and can be theoretically described within statistical 
mechanics. Given the interaction between the species (particles, spins, etc.) 
the partition function can be calculated, in principle. Its logarithm yields, e.g. 
for a canonical ensemble, the free energy which is singular at the equilibrium 
transition point. This allows to fix this point from first principles. 

Besides such order-disorder transitions there also exist transitions between 
disordered phases. Excluding liquid-liquid transitions these will be called glass 
transitions. A prominent example is the transition from a supercooled melt 
of Si02-molecules to an amorphous phase which is the well-known "window 
glass". One distinguishes two types of glass transitions: spin glass transitions 
PI E] and structural glass transitions [HI IH 13 EJ- Their main difference is 
that the former occur mostly in systems with quenched disorder and for the 
latter the disorder is self-generated. We stress that this classification should 
not be taken too strict since there are also models without quenched disorder 
showing spin glass behavior [7j. Typical systems undergoing a structural glass 
transition are liquids, particularly molecular liquids, like the famous example 
of SiC>2. In recent years investigations of specific spin glass models, e.g. Potts 
glass and so-called p-spin models with p > 3 where p spins are coupled by 
randomly frozen-in (i.e. quenched) infinite range interactions, have revealed 
some similarities with structural glasses [E|. We will come back to this point 
in the 4th chapter. 

In contrast to conventional order-disorder transitions, glass transitions are 
less well understood. There is a broad consensus that a spin glass transition 
exists and that an appropriately disorder-averaged free energy is singular at 
the transition point in case of mean field models, i.e. models with infinite 
range interactions or infinite dimensions. But this is less obvious for short 
range interactions. The situation for the structural glass transition is even 
less satisfactory, although substantial progress has been made in the last two 
decades. 
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In this article we will mainly focus on the structural glass transition. Spin 
glass behavior and spin glass transition will be discussed in this monograph 
by H. Horner. For further details and references on spin glasses, the reader 
may consult his contribution. 

As mentioned above, structural glasses can be obtained by cooling a liq- 
uid. In order to bypass crystallization one has to choose a finite cooling rate. 
For good glassformers, like SiC-2, this rate can be rather modest whereas bad 
glassformers, like most metallic glasses, require extremely high cooling rates. 
Under such a cooling process the shear viscosity rj(T) increases. Close to the 
so-called calorimetric glass transition temperature T g the supercooled liquid 
falls out of equilibrium and becomes a glass. At T g thermodynamical quan- 
tities like density n(T), specific heat c p {T) at constant pressure, etc. show 
cross-over behavior, i.e. the slope of n{T) and c p (T) make a more or less 
well pronounced jump, depending on the cooling rate. T g itself depends on 
the cooling rate, too. Although T g plays an important practical role, it is 
less interesting from a fundamental point of view, due to its cooling rate 
dependence. Besides T g there are at least three more characteristic temper- 
atures T , Tk and T c . For many glass formers, rj{T) can be fitted by the 
Vogel-Fulcher-Tammann law 

r,(T) = Cexp A , T > T (1) 

k B (T - lb) 

with A > and C > 0. The shear viscosity diverges at To. Extrapolat- 
ing the excess entropy S excess (T) of the supercooled liquid with respect to 
the crystalline phase to lower temperatures there is the so-called Kauzmann 
temperature Tk at which 5 excess vanishes: 

•SexcessfTA') = . (2) 

Since it is argued that a disordered phase should not have a smaller entropy 
than the crystalline one, S excess can not become negative. Therefore, the sys- 
tem has to undergo a static glass transition at Tk- This conclusion, however, 
is not compelling, since there exists inverse melting, i.e. liquids freeze when 
heated or crystals melt when cooled [2j . In that case the total entropy of the 
crystal is higher than that of the liquid. A recent discussion of the Kauzmann 
problem can be found in Ref. ^H]- T 9 ,Tq and Tk have played an essential 
role for many decades. In 1984 quite a new theoretical approach, the mode 
coupling theory IH\> has shown that there is a critical temperature T c , at 
which a dynamical glass transition takes place. One of the main features is 
that the nonergodicity parameters f(q,T) which can be considered as glass 
order parameters change discontinuously at T c : 

/(«.*)={> 2 ; t<tI . ( 3 ) 

Since then numerous experimental investigations and computer simulations 
were stimulated (see reviews ^2 E] and Ref. ^31 El)- They have shown 
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new characteristic dynamical features close to the dynamical glass transition 
point T c consistent with mode coupling theory (see also the contribution by 
U. Buchenau in this monograph). 

This short exposition of some of the characteristics of glassy behavior 
should have given a first impression on how diverse the phenomena in the 
glass transition region can be. Therefore it is obvious that a successful theo- 
retical description which covers all facets is extremely hard. There are mainly 
two possible theoretical approaches: phenomenological or microscopic ones. 
Phenomenological theories start from some of the phenomena of glasses, and 
are named thereafter. Based on these phenomena a theoretical description 
is developed capable of describing the observed phenomena. In several cases 
appealing "physical pictures" are used. However, a couple of assumptions are 
made which are not proven. The predictive power of such phenomenological 
approaches is rather limited. This is quite different from a microscopic theory. 
By microscopic we mean that the physical quantities can be calculated from 
first principles if the interactions between the species are given. Since the 
glass transition region is located at rather high temperature quantum effects 
can be neglected. Therefore, a microscopic theory starts from a classical TV- 
body problem. The next chapter will discuss some of the phenomenological 
theories. The major part of this article is devoted to microscopic theories: 
The 3rd chapter describes mode coupling theory and the 4th chapter the 
replica theory for structural glasses. 

Finally we want to stress that the present contribution presents a selection 
and does not aim to be complete. This holds mainly for the phenomenologi- 
cal models. We also do not discuss the potential energy landscape"potential 
energy landscape" approach W\ which recently has led to new interesting re- 
sults |17II18| . Almost all of them were obtained from computer simulations. It 
would be desirable to complement these investigations by analytical theories. 

2 PHENOMENOLOGICAL APPROACHES 

The presentation in this chapter will be rather short. More details and addi- 
tional phenomenological approaches can be found in the monographs El El 
and in the review 0]. 

2.1 Adam-Gibbs Theory 

In 1965 Adam and Gibbs suggested a theory based on the assumption that 
dynamically cooperative regions occur when decreasing the temperature to- 
wards the glass transition point The particles in these regions perform 
cooperative motion, which leads to a reduction of the configurational degrees 
of freedom. Here we follow the description of Adam-Gibbs theory as given in 
Ref. 0J. The basis of that theory is a number of assumptions: 
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1) For given temperature T there are dynamically cooperative regions la- 
beled by 1, ... ,i, ... , n(T) with N±(T), . . . , A^(T), . . . , N n ^(T) particles. 
The number n(T) of regions decreases with decreasing temperature (see 
Fig.0. Due to the conservation of N, the total number of particles, it is: 

n(T) 

^7V 4 (T) = iV . (4) 

1=1 

2) These regions take two configurations, only. Accordingly the entropy s 
per region is: 

s = fc B ln2 . (5) 

The assumptions that each region takes two configurations, only, is not 
crucial, but it should be a finite number of order one. 

3) Fluctuations of Ni(T) are small, i.e. it is: 

N t (T) w iVo(T) = N/n(T) . (6) 

because of Eq. Q. 

With these three assumptions we can relate the average number Nq(T) of 
particles within a dynamically cooperative region to the configurational 
entropy: 

S C (T) — ks In [number of total configurations] . (7) 

With 2™( T ), the number of total configurations we get from Eq. by use 
of Eq. ©: 



THEORIES OF THE STRUCTURAL GLASS TRANSITION 5 

N Q (T) =Nk B \n2/S c (T) (8) 

i.e. Nq(T) and therefore the size of those regions is inversely proportional 
to the configurational entropy. This is plausible since the number of pos- 
sible configurations decreases with increasing Nq(T). 
The next important steps are to assume: 
4) There is a temperature Tk at which Nq(T) becomes infinite (for N = oo), 
Then Eq. JHJ implies 



Sc(T) 



> 
= 



T>T K 
T<T K 



(9) 



i.e Tk is the Kauzmann temperature |2fl] . 
5) The transition between both configurations of a region is an activated 
process with an activation energy: 



E(T) = e N (T) 



(10) 



where e may be weakly T-dependent. Then the transition time t(T) is given 
by: 



r(T) ~ eMPE(T)) = exp(/3 eo 7V (T)) 
Substituting Nq from Eq. JSJ, we arrive at: 

"e iVln2" 



t(T) ~ exp 



(11) 



(12) 



TS C {T) 

Assuming that S(T) vanishes linearly at Tk we obtain from Eq. J12II the 
Vogel-Fulcher-Tammann law close to Tk- 



t(T) ~ exp 



e iVln2 



(13) 



T K S' C {T K ){T -T K )_ 
with To = Tk- 

This example demonstrates that the appealing picture of dynamically coop- 
erative regions in combination with several assumptions allows to derive the 
Vogel-Fulcher-Tammann law. But it also shows that no additional predictions 
are made. Furthermore, it is obvious that the growth of the dynamically co- 
operative regions would be accomplished by a divergent length scale, which, 
however, has never been found in experiments or computer simulations. 



2.2 Free-volume theory 

This phenomenological description has been made by Cohen and Turnbull 
in 1959 I n order to illustrate their idea we choose a system of hard 
spheres with average density n. Let us fix the positions of all spheres except 
of the i-th sphere. Then the z-th sphere can move freely in the so-called free 





Fig. 2. Illustration of the free volume of one sphere (full circle) within a fixed 
configurations of other spheres (open circles). The free volume is shown as hatched 
regions. Upper panel: T > To for which the free volume forms a percolation cluster. 
Lower panel: T < To for which the free volume does not percolate but gets localized. 

volume Vf(i) (see Fig. [2J. Although Vf(j) is correlated with v/(i) for i ^ j 
one assumes: 

1) Every sphere has a free volume Vf. 

i -> v f (i) (14) 

2) Vf(i) are independent random numbers with probability density p(vf) 
which is assumed to be exponential: 



p(vf) = A/"exp(— Vf/vf) 



(15) 
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with the mean free volume per particle: 

1 N 

i=l 

where N is a normalization constant. 
3) It is obvious that Vf(i) and therefore Vf decreases with increasing den- 
sity or decreasing temperature (in case of "soft" particles) . An essential 
assumption is that there is a temperature T > such that the mean free 
volume vanishes at T : 



u A T )-{ , T<T Q W 

with a the expansion coefficient. 
4) The inverse shear viscosity is proportional to the probability that the free 
volume per particle is larger than a certain value , i.e.: 



V l i T )~ / dv fP( v f) 



(18) 



Substituting Eq. {T^J an d Eq. {T^Jl into Eq. {Hi we get the Vogel-Fulcher- 
Tammann law: 



r)(T) - exp 



J f 



a{T-T )_ 



(19) 



2.3 Extended free-volume theory 

The free volume theory has been extended by Cohen and Grest [22 . The 
main idea of this extension is to make a connection between the glass transi- 
tion and percolation. If the density n (temperature T) is below n g (above T g ) 
the individual free volumes Vf(i) overlap or, in other words, the free volume 
percolates in such a way that a (macroscopic) percolation cluster exists (see 
upper panel of Fig.EJ. In this case a particle can move macroscopic distances, 
i.e the diffusion constant is finite. Now, increasing n (decreasing T) there may 
exist a critical value n g (or T g ) at which the percolation cluster disappears 
which implies that the particles become localized in a glass phase with zero 
diffusion (see lower panel of Fig.[2j|. If this scenario would be correct the glass 
transition point would coincide with the percolation threshold of the free vol- 
umes Vf(i). Since percolation is well understood this relationship would 
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imply several characteristic features, e.g. the fractal nature of the percolation 
cluster, close to n g (or T g ). However, there is no experimental evidence for 
such a fractal behavior. The main critique on the free volume theory is that 
the density used in the lower panel of Fig. [2 to demonstrate the localization 
of the free volume still corresponds to a liquid. Computer simulations show 
that the density of a liquid of e.g. hard discs is already higher than that used 
in the lower panel of Fig. [21 Accordingly, the free volume, as defined above, 
is already localized in the liquid phase. 

2.4 Gibbs-DiMarzio theory 

This approach is not really phenomenological. But the model which has been 
studied by Gibbs and DiMarzio is very special and has a limited range of ap- 
plicability. Therefore it is included in this chapter. The model by Gibbs and 
DiMarzio is capable to explain the vanishing of the configurational entropy 
|24] . A system with N polymers on a cubic lattice with N s lattice sites is con- 
sidered. Each polymer consists of I monomers and can form conformations 
labeled by 1, 2, ■ • • . Let n, be the number of polymers with conforma- 

tion i and P({m}, N, N s ) the probability that N polymers on a lattice with 
N s sites built a set {rii} of conformations. Associating an energy e% with each 
individual conformation the total energy E is given by: 

E = Y / n l e l . (20) 

i 

Then the number of configurations f2(E,N,N s ) with energy E is: 

Q(E, N, N s ) = P{{mh N, N s ) 5{E - £ sm) (21) 

w 

from which one obtains the configurational entropy: 

S C (E,N,N S ) = k B lnf2(E,N,N s ) . (22) 

Of course, the nontrivial problem is to determine P and to perform the sum 
over {m} in Eq. (J2lj. Making use of a mean field approximation due to Flory 
and Huggins j2^] one finally gets a critical energy Ek, which corresponds to 
the Kauzmann temperature Tk with: 

; lit . <m 

Since mean field approximations tend to produce phase transitions even for 
systems which do not show such transitions, the validity of results Eq. ffify 
is not obvious. 
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3 MICROSCOPIC THEORY: MODE COUPLING 
THEORY 

We hope that the short presentation in the 2nd chapter has demonstrated 
that the status of those phenomenological theories is not satisfactory. Al- 
though, plausible "physical pictures" are involved they are based on a couple 
of crucial assumptions. However, the validity of those remain completely un- 
clear. This demands for a microscopic approach based on first principles. 
For the first time, such an approach to structural glass transitions was made 
in 1984 by Bengtzelius, Gotze and Sjolander [H]. Starting from the liquid 
side these authors applied the mode coupling theory (MCT), developed by 
Kawasaki [26 in order to describe the critical slowing down close to a critical 
point, to the relaxation of the density fluctuations in a supercooled liquid. 
The result is an equation of motion (see below) for the intermediate scatter- 
ing function S(q,t) for a simple liquid which can be measured by neutron- 
and light scattering (see the contributions by U. Buchenau and J. B. Suck in 
this monograph). The properties of the solution S(q, t) of the MCT-equations 
were mainly investigated in great detail by Gotze and his coworkers. These 
results can be found in the reviews [23 EEl EE El Ei| and in their references. 
Since most glass formers are molecular systems which involve translational 
and rotational degrees of freedom, MCT had been extended to a single linear 
molecule in a simple liquid by Franosch et al. [201 and to a liquid of linear 
molecules and arbitrary molecules by Scheidsteger and the author jHI] and by 
Fabbian et al. JUL respectively. This was accomplished by the use of the ten- 
sorial formalism which allows the separation of translational and rotational 
degrees of freedom. Alternatively, one can also use a site-site representation 
for molecular systems. MCT based on such a description was worked out by 
Chong and Hirata [22] and Chong, Gotze and Singh |33j . 

Before discussing MCT let us anticipate that the derivation of the MCT- 
equations requires some more or less strong approximations. Although, these 
approximations can not be controlled, e.g. due to the lack of a smallness 
parameter, it is interesting to notice that the MCT-equations for a simple 
liquid were also obtained by quite different approaches which are the use 
of generalized fluctuating hydrodynamics by Kirkpatrick 84j and Das and 
Mazenko [35^ , density functional theory by Kirkpartick and Wolynes |32] and 
recently by the use of the equation of motion for the microscopic density 
p{q, t) in conjunction with assuming the density fluctuations to be Gaussian 
by Zaccarelli et al. [2Zj. That these quite different approaches lead to the 
same mathematical structure of the MCT-equations shows its "robustness". 
In addition, it has been argued that the MCT-equations become even exact 
in infinite dimensions [HEj- Further support for MCT comes from a spherical 
model (without quenched disorder) [38.; and for spin glass models with infi- 
nite range interactions (see Ref. [B, .'59 and their references) for which the 
corresponding q- independent MCT-equation has been proven to be exact. 
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Now we will turn to the discussion of MCT. This will be done in two 
sections. In the first one we present the derivation of the MCT-equations and 
in the second one the properties of their solutions. 

3.1 Derivation of the MCT-equations 

We will restrict ourself to a simple liquid of N identical particles with mass 
m in a finite volume V. We assume two-body interactions v(x — x') such that 
the classical hamiltonian is given by: 

H({ Xi }, { Pi }) = ]T ^vl + V({x t }) (24) 

n 

with the potential energy: 

V({x i }) = ^v(x i -x j ) ( 25 ) 

{xi} and {pi} are the corresponding positions and momenta, respectively. 
Fixing an initial point {x^, {p^ in phase space, the phase space point 
{x{(t)}, {pi(t)} at time t is determined by Newtonian's equation of motion 
or equivalently by the hermitean Liouville operator: 

Xi{t) = e iCt Xi , Pi (t) = e iCt Pi . (26) 

Of course, one can not solve these equations of motion for a macroscopic 
system. Therefore one has to use a theoretical framework which restricts 
itself on the relevant variables. For a liquid this is the microscopic density: 

p{x,t) = ^6{x-x n {t)) (27) 

n 

where the dependence of the positions x n (t) on the initial point {x{\, {p{\ 
is suppressed. Now we recall that p(x,t) will have a slowly varying part in 
the strongly supercooled regime, besides fast motions (vibrations) around the 
quasi-equilibrium positions. Therefore, we can consider p(x,t) or its Fourier 
transform 

p(q,t) = Y,e iqXnit) =e tCt p(q) (28) 

n 

as a slow variable. If p(qr, t) is slow, then the current density 

n 

is slow, too, since it is related to p(q, t) through the continuity equation: 
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p(q,t) = i£p(q,t) = iq-j(q,t) . (29) 
We can continue taking time derivatives. Then p(q, t) is given by: 

p(q, t)=iq-^T— p n (t) e iqx " w + kinetic part . (30) 

n 

Using p n /m = —dV({xi})/dx n and that Eq. can be 
rewritten as 

V({xi}) = i Jd 3 x J d 3 x' v(x - x') p(x)p(x') 

= ^E^(")^) ( 31 ) 

q 

it is easy to prove that: 

p{Q, t) = 5Z'^o(q, 9i, qi) P(qi,t) p(q 2 , t) + kinetic part (32) 
with the bare and time-independent vertex: 



^0(9,91,92)= q ■ [qxv{qi) + 92-0(92)] • (33) 
Here we have used that the Fourier transform v(q) of the pair potential 



is real and denotes summation such that qi+q 2 = 9- The result Eq. l|32ll 
reveals that the "force" p(q, t) contains contributions from a pair of modes, 
which are coupled by the bare vertex wq. It is rather obvious that the use of 
a r-body interaction would yield a contribution p(qi,t)p(q2, t) ■ . . . ■ p(q r , t). 
Consequently, the "force" p(q, t) must be slow, as well. By continuing this 
procedure one obtains a set of slow variables. In the following, however, we 
will restrict ourself onto the two variables p(q, t) and j(q, t) = q j(q, t)/q(q = 
\q\), the longitudinal current density. But we have to keep in mind that 
p{q, t) also is slow. Having chosen p(q) and j(q) (at t = 0) as slow variables 
one can apply the Mori-Zwanzig projection formalism [401 l4"T] to derive an 
exact equation of motion for the normalized intermediate scattering function 
$(q,t)=S(q,t)/S(q) with 

S(q,t) = ±(p(q,typ(q)) = ±(p*( q )e~* ct p{q)) (34) 

where we used the hermiticity of C. S(q) is the static structure factor which 
depends on the thermodynamical variables T, n, etc. The exact Mori-Zwanzig 
equation for <P(q,t) reads: 

t 

S(q, t) + Q 2 q *(g, t) + J d£ M(q, t - t')4>{q, t') = (35) 
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M(q,t) = - — -(ptqyQe-^Qpiq)) (36) 
and the microscopic frequencies: 

W^^)" 2 . (37) 

The initial condition is <P(q, 0) = 1, and ^(<7,0) = for all q. The result 
Eq. I|35ll makes obvious that the problem to calculate <P(q, t) has been shifted 
to the calculation of M(q,t), which seems to be hopeless, as well. But this 
is not really true. In contrast to $ it is possible to approximate the memory 
kermel M . M(q, t) is the correlation function of the "forces" Q p(q), however, 
with the reduced Liouvillian £' — QCQ. Q is the projector (see below) 
which projects perpendicular to both slow variables p{q) and j(q). Since 
p{q) contains a coupled pair of modes p{qi) p{q2) (cf. Eq. i!32|) ) and because 
Qp(qi) p{q2) 7^ 0, the "force" Qp{q) still contains a slow part. This suggests 
to use the following approximation |27] : 

Qp(q) -h. Q\p(q)) » VQ\p{q)) (38) 
where V is the projector onto pairs of modes: 

V= 9(Qi^;qW2)\p(qi)p(q2))(p(q' i rp(q , 2r\ (39) 

HQ2,q' 1 q' 2 

and g is determined such that V 2 — V. The reader should note that we have 
introduced a bra- and ket-notation (| and |), respectively, like in quantum 
mechanics. This can be done since the canonical average (A*B) of two phase 
space functions A* and B can be interpreted as scalar product (A*\B). It is 
the existence of this scalar product which allows to introduce projectors. Sub- 
stituting Eq. <EBl and Eq. (jSHJ into Eq. ll37)j) relates M(q, t) to the correlation 
function 

(p(q' 3 rp(q' 4 y\e- iC ' t \p(q 1 )p(q 2 )) . (40) 

This relationship involves g(qiq 2 ; q^q^gilzl^ QWd and (p(q)*Q\p(q3)p(q4)} 
(p(q' 1 )* p(q' 2 )*\Qp(qi)} which are static quantities. Now, the crucial approx- 
imation is the factorization of the correlator Eq. ll4Ti)l and simultaneously 
replacing £ by C: 



(p(q' 3 Tp(q' 4 r\e- iC ' t \p(qi)p(q2)) » 
« [(p( qi r\e- iCt \p( qi )) (p(q 2 y\e- iCt \p(q 2 ))5 q , qi S q , q2 + 

+ (1 — 2)] 

= N 2 S(q i )S(q 2 )^(q 1 ,t)^q 2 ,t){S q/3qi S qiq2 +S qkq J qtiqi ] . (41) 
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The condition V 2 = V implies that g(qiq 2 ', Q1Q2) ls the inverse of the 
"matrix" ((p(q'i)* p(q' 2 )* p(<li) pilz))) ■ Since this "matrix" equals the correlator 
Eq. ijlflj) at t = it is approximated by the r.h.s. of Eq. (fill) , at t = 0. Using 
this approximation one immediately finds: 

qi 3q' 2 q2 + ^<3i<j2 <W ■ (42) 

UsingQ=l-7' = l-(iV5( 9 ))- 1 |p( 9 ))<p(g)*|-(iVfc B r/TO)- 1 |j(g))(j(g)1 
the static correlation (p{q)* Q\p{qi)p{q 2 )) can be expressed as follows [27| : 

(p(g)*0|p(gi)p(«2)) - 5(gi) 5(g 2 ){-ng • [ qi c{q x ) + g 2 c(g 2 )] + 

m 

+ q 2 [l-(p(qy P {q 1 )p{q 2 ))/{NS{q)S{q 1 )S{q 2 ))]}5 qi+q2 , q (43) 

where we introduced the direct correlation function c(q) defined by S(q) 
I - nc(qf)] -1 . Substituting Eqs. QEJ-O into Eq. ® yields finally the 
MCT-approximation for the memory kernel: 

M(q, t) « 1/, tf(t) + /2 2 m(g, t) (44) 

= ^7 ^'^(g,gi,g 2 )^(qi,i)<?(g2,i) (45) 

qi,q2 

with the positive vertices: 

V(q, qi,q 2 ) = - 5(g) 5(gi) 5(g 2 ) q~ A {nq ■ [gic(gi) + g 2 c(g 2 )] - 
n 

-9 2 [1 - {p{q)*p{qi)p{q 2 ))/{N 5(g) 5(g x ) 5(g 2 ))]} 2 . (46) 

The first term on the r.h.s of Eq. 144t accounts for the fast part of the 
"force" Qp{q) leading to a frictional contribution. 

The equations Eqs. jSHt i CUt and lH4l - (Hfil) are the mode coupling 
equations. Due to the MCT-approximation they are a closed set of integro- 
differential equations for the normalized correlator <P(q,t) with initial con- 
ditions <P(q,0) = 1 and ^(g,0) = 0, because of time reversal symmetry. 
As an input they only need the static two-point correlator 5(g) (or equiva- 
lent^ the direct correlation function c(g)), the static three point correlator 
(p(q)*P{Qi) p{q2)) and the frictional constants v q . v q which can only be deter- 
mined from kinetic theory does not influence the glassy behavior. Therefore 
it can be put to zero. The remaining static two- and three point correlators 
5(g) = {p{q)* p{q)) / N and (p(q)* p(qi) p[q%)) can be calculated for given po- 
tential energy V({xi}). This is what makes MCT a microscopic first principle 
theory. Here a comment is in order: MCT will be applied to the supercooled 
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liquid, i.e. to a temperature regime where the stable thermodynamical phase 
is a crystal. In order to study the glass transition one has to use the static 
correlators for the supercooled liquid and not for the crystal. MCT allows to 
predict the time dependence of the density fluctuations provided both static 
correlators are known. They can be obtained either from analytical approxi- 
mation schemes 0Hj or from experiments and simulations. Application of the 
convolution approximation |4fl| : 

(p(g)*p(<?l)p(g 2 )) « N S(q)S( qi )S(q 2 ) S qi+q2 , q (47) 
leads to a further simplification of the vertices Eq. ll4(ill : 



V(q, qi,q a ) = nS(q)S( qi )S(q 2 ) [q ■ (qic(qi) + q 2 c(q 2 ))] 2 /< Z 4 (48) 

which involve S(q) (or c(q)), only. For SiC>2-liquids it has been demonstrated 
by Sciortino and Kob ^2 that a satisfactory agreement of, e.g. the critical 
nonergodicity parameters determined from a MD-simulation with those from 
MCT is only obtained with the vertices from Eq. I|43|) where the three-point 
correlator is not factorized. This is rather plausible, because SiC-2 is a cova- 
lent glass former with bond orientational correlations which are completely 
neglected by the convolution approximation Eq. I|47|) , 

3.2 Solutions and predictions of MCT 

The main question which arises is: How can one detect a glass transition 
within MCT? The answer is rather simple. Let us use the nonergodicity pa- 
rameters defined by: 

f(q) = lim $(q,t) = - lim zS(q,z) . (49) 

t — >oo z — too 

These parameters, which are just the infinite time limit of <£(q, t) or the zero 
frequency limit of its Laplace transform 

oo 

S(q,z) = i J dt<P(q,t)e lzt , Im z > (50) 
o 

can be used as glass order parameters, because they vanish in an ergodic phase 
(= liquid phase) and are nonzero in a nonergodic one, which is interpreted as 
a glass phase. To be more precise, this is true for the correlation function of 
the density fluctuations Sp(q,t) = p(q,t) — (p(q,t)), only. Since 5p(0,t) = 
and Sp(q,t) = p{q,t) for all q ^ 0, it is sufficient to investigate <P(q,t). Thus 
it is: 

/«"={>» : < 51 » 
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Since <P(q, t) — > for t — > oo, it is easy to show that Eq. j^H yields a nonlinear 
set of algebraic equations for f(q): 

Jiq] =nf(q)} (52) 



1-/(9) 
with: 

nf(9)] = 27 E'^9i,g 2 )/(9i)/(9 2 ) . (53) 

9i,9a 

One can prove that the long time limit f(q) is distinguished from other 
possible solutions of Eq. <l52l) , (Ifi3l) by the following properties [22j 

(i) f(q) is real (since $(q,t) is real) 

(ii) < f(q) < 1 

(iii) If several solutions fi(q) > fziq) > h{q) > •• ■ exist, the long time limit 
is the largest one: f(q) = fi(q) ■ 

The positiveness of the vertices is crucial for (iii) . 

It is obvious from Eqs. JH2J, that f(q) = is always a solution. 
In order to check whether a nontrivial solution exists, let us neglect the q- 
dependence for a moment. This leads to a so-called schematic model the 
.TVmodel, where !F[f] = vf 2 , v > 0, i.e. we have to solve 

7^7 =«/ 2 • (54) 
For / ^ this leads to a quadratic equation with solution: 

h/2 = \[l± Vl-4/«] . (55) 

Since fx/2 is n °t real for v < v c = 4, the only physical solution in this 
range is /o = 0. For v > v c two real solutions /i > fa bifurcate from the trivial 
one. This scenario is illustrated in Fig. El Hence, we get for the i<2-model 

/(v) = U[i + v/T^H , v>v c ^ 

i.e. at the critical control parameter v c = 4 the nonergodicity parameter 
changes dis continuously from zero to the critical nonergodicity parameter 
/ c = f(v c ) = 1/2. This type of transition is called type-B transition, in con- 
trast to a type- A transition which is continuous [27]. Such a type- A transition 
occurs, e.g. for the .Fi-model: F\f] = vf, v > 0. 

This simple example has shown that an ergodic-to-nonergodic transition 
can occur at a critical coupling constant v c . Since the vertices V(q, qi, q?) are 
positive and < f(q) < 1, !F[f(q)] will diverge in the strong coupling limit 
^(<7j <7ij 92) — * 00. Therefore Eq. (lo2fl can only be fulfilled if the nontrivial 
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f(u) 
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fo(u) 



^(u) 




u 



Fig. 3. Qualitative control parameter dependence of the nonergodicity parameter 
of the J-2-model. /o(i>) = is the trivial solution. At v c = 4 two new solutions fi(v) 
and /2(f) bifurcate. 

solution converges to one. Therefore there must be a critical hypersurface in 
the control parameter space at which transitions from /(g) = to /(g) ^ 
happen. Because the vertices depend through S(q) on the thermodynamic 
variables T, n etc. and increase with decreasing T (increasing n) there will 
be a critical temperature T c (critical density n c ) at which the system under- 
goes an ergodic-to-nonergodic transition, i.e. a glass transition. From this we 
can conclude that MCT yields a dynamical glass transition whereas static 
quantities, e.g. S(q), are not singular at T c (or n c ). Having established the 
existence of a glass transition singularity, one can study the dynamics close 
to it. For large times (small frequencies) compared to the microscopic time 
scale I2~ 1 (i7 g ) the Laplace transform of Eq. (l35l) yields with Eq. lEjl : 



Here we used v q = 0. The reader should notice that Eq. 1)5 7JI does not in- 
volve anymore the microscopic frequencies Q q . The schematic representation 
of a solution of the MCT-equations Eqs. l(35)l . l(44ll - lllfill . including the mi- 
croscopic time regime, is shown in Figure^l This Figure clearly demonstrates 
the existence of a critical temperature T c at which a jump of /(g) occurs. 
The generic behavior close to T c is given by (cf. also Eq. ®) : 




, T > T c 

f c {q) + const(T c - Tf/ 2 ,T<T C 



(58) 
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«Kq,t) 



1 




f c (q) 



= 0.01 



In t 



Fig. 4. Qualitative time- and temperature-dependence of the intermediate scat- 
tering function $(q,t) for fixed q. The values 0.1 and 0.01 for (T - T c )/T c should 
not be taken literally. They should only indicate that T is chosen closer and closer 
to T c . 

where the constant is g-dependent. We can observe from FigureQJthat close 
to T c there is a time scale on which the correlator <P(q, t) is close to the critical 
nonergodicity parameter f c (q) = f(q)\ T=Tc - This suggests the Ansatz: 



with \G{t)\ -c 1. It is important to realize that the q- and ^-dependence is 
factorized. This is related to the type of bifurcation scenario for which the 
non-degenerated largest eigenvalue of the stability matrix of the linearized 
Eq. llo2ll becomes one. Therefore at T c only one unstable eigenvector occurs. 
The critical amplitude h(q) is the amplitude of that eigenvector. Substituting 
Eq. 15911 into Eq. Q57JI and expanding up to quadratic order in G one obtains: 



$(q,t) = f c (q) + h(q)G(t) 



(59) 



(7 + \{-zLT[G 2 {t)]{z)} 



{-zG(z)} 2 = 



(60) 



with the separation parameter 




(61) 



(ctq is positive) and the exponent parameter 




(62) 



The separation parameter is a measure of the distance from the transition 
point. The importance of A will become clear below. Explicit expression for 
a and A are given in Ref. |27j . 



18 Rolf Schilling 



At the transition point, i.e. for a = 0, the solution of Eq. (j6fi)) is the 
critical law: 

G{t) ~ r a . (63) 
where the exponent a is a solution of: 

r 2 (l-«) 

with r{x) the Gamma function. Since A is uniquely determined by S(q) at 
T c it can be calculated microscopically, from which a follows, a is restricted 
to0<a<l/2. For a = a higher order bifurcation scenario occurs |27] , 
To solve Eq. 116011 for T ^ T c it is convenient to introduce a correlation scale 
V^by: 

G(t) = ^\ g ± (t/t a ) , a^O (65) 

where the time scale t a and the master functions g± can be obtained as 
follows. Introducing Eq. Illiot into Eq. KiOII and using t = tjt a and z = t a z 
leads to: 

±1/1 + XLT[gi(t)](z) + z(g±(z)) 2 = (66) 
which has the solution: 

S±(*)~r» (67) 

for t <C 1, i.e. t <C t a , since ±1/1 can be neglected due to z 3> 1. Substitution 
of Eq. ll6T|) into Eq. l(65|l and taking into account that G(t) must reduce for 
<7 — > to the critical correlator Eq. Ki.'ill allows to fix t a ; 

t a (T)~\ a \-h ~\T-T C \-^ , T^T C . 
For i > 1 one obtains E3: 



with B > 0. 

Again, the positive exponent b is determined by the exponent parameter 

A: 

_ r\i + b) 
x r(i + 2b) • (70) 
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t <Kq»t) 




Fig. 5. (S>(q,t) for fixed q and a single temperature T > T c (solid line). Both 
power laws, the critical law and the von Schweidler law, are shown as dashed lines. 
The horizontal dashed line is the plateau height equal to the critical nonergodicity 
parameter / c (<?) 



In the liquid phase (minus sign in Eq. {HHJl) we obtain, besides the critical 
law Eq. Ijfi.SJI a second power law, the so-called von Schweidler law. Substi- 
tuting g~(t) from Eq. ijfifljl into Eq. leads to 

G(t) ~ (t/r(T)) b (71) 

with the second time scale: 

T (T)~\*\-&t a ~(T-T c )-'', T>T C (72) 

and the exponent 

i = k + k ■ <w 

The von Schweidler law, Eq. {ZD, is van d for i > 1, i.e. t > t a {T) and 
t <C r(T) because of \G(t)\ <§; 1. Both power laws are shown in Figure El 
The critical law describes the relaxation to the plateau value f c (q) and f(q) 
above and below T c , respectively, and the von Schweidler law the initial decay 
from the plateau for T > T c , Both time scales t a and r exhibit power law 
divergence at T c . t diverges faster than t a , due to 7 > /l(2a) (see Fig- EJ) - 

These fascinating results proven by Gotze in 1984 05 an d 1985 jT4 were 
absolutely new. The time range to <C t -C r(T) (to is a microscopic time 
scale ~ ^q 1 ) which includes both power laws have been called (fast) (5- 
relaxation. It has a simple physical explanation. If the temperature is low 
enough each particle feels a cage on a time scale t a . The first relaxation 
step (critical law) is a relaxation within the cage and the second one (von 
Schweidler law) is related to the "opening" of the cage for t > which 
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relaxation time 



t, 




T 



Fig. 6. Qualitative T-dependence of both relaxation times t a (T) and t(T) which 
exhibit power law divergence at T c 

is the initial stage of the structural relaxation. The fact that the q- and t- 
dependence of the correlators factorize in the /?-regime is one of the highly 
nontrivial predictions of MCT which has not been found before in condensed 
matter physics. Numerous experiments ( see a l so the contribution by U. 
Buchenau in this monograph) , and computer simulations jl.'U 1141 IT5] have 
given consistent results with regard to these MCT-predictions. 

What remains is to study the dynamics for t of the order or much larger 
than t(T) which is called a-relaxation. There the factorization does not hold 
anymore. Therefore one has to solve Eq. (JSJJ. This can only be done nu- 
merically. But there is one important feature, which is the scale invariance of 
Eq. I|57jl. This means that a scale transformation t — > t ■ y, <P(q, t) — > <P(q,t-y) 
or z — > z ■ y~ x , <P(q, z) — > y^(g, z ■ y^ 1 ) leaves Eq. if^Tjl invariant. This allows 
to introduce a g-dependent master function & q (i) such that 



i.e. the time and temperature dependence of <P appears only in the com- 
bination t/r(T). Eq. lT74"l) is well-known in the glass science as the time- 
temperature superposition principle. The master function <£ g (i) can be de- 
termined from a numerical solution of Eq. ll57l) . The validity of Eq. ij^Tjl and 
Eq. I|74H has also been test experimentally an d by computer simulations 



*(g,t;T) = «P,(t/r(T)) 



(74) 



Gain]. 
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4 MICROSCOPIC THEORY: REPLICA THEORY 

In 1996 Mezard and Parisi H5j have presented a replica theory for the struc- 
tural glass transition. This theory has been inspired by spinglass theory P3E]. 
Kirkpatrick, Thirumalai and Wolynes first noticed the analogy between sp- 
inglass and structural glass transitions In 1987 they have shown that 
mean-field spin glass models with a one-step replica symmetry breaking j2] 
exhibit a dynamical (MCT-like)glass transition at T c and a static one at T s , 
which is below T c . At T 8 the configurational entropy vanishes. Hence T s can 
be identified with Tr-, the Kauzmann temperature. In the next section we 
will describe the physical picture behind this theory and its thermodynamic 
formulation. The second section then contains the microscopic theory which 
can be considered as a first principle approach of an earlier attempt by Singh, 
Stoessel and Wolynes to describe the glass transition by a density functional 
theory |47i| . 

4.1 Thermo dynamical description 

Let p{x) be the local particle density of a liquid with TV identical particles 
in a volume V . A theorem 02] guarantees that there exists a free energy 
functional J-[p(x),T] such that the equilibrium phases are obtained from 

Sp{x) 

J-[p(x),T] is not known exactly. In many cases one uses an approximation 
suggested by Ramakrishnan and Youssouff [4Hj. Let p^ a \x) (T-dependence is 
suppressed), a = 1, 2, 3, ... be solutions of Eq. lT75l) (local minima!) with free 
energy per particle f a (T) = T[p( a \T]/N and let A/(/, N, T) be the number 
of solutions with free energy / = T '/N at T. The configurational entropy per 
particle S c (f,T) (cf. 1st chapter) is defined by: 

Af(f,N,T)=e X p[NS c (f,T)] . (76) 

At high temperature the equilibrium phase is given by the uniform density 
solution p{x) = n = N/V of Eq. lf75l) . Inspired by mean field spin glasses 
with a discontinuous transition, Mezard and Parisi 001 HO] assume that at 
the MCT-temperature T c an exponential number of solutions p^ a \x) occur 
for / between f mi „(T) and / max (T), i.e. 

S o (/,T)>0, / min (T) < / < / max (T) . (77) 

Above f ma ,x{T) and below f m i n (T) it is S c (f,T) — 0. This situation is il- 
lustrated in Figure S c (f,T) varies smoothly with T and is concave in /, 
i.e. d 2 S c {f,T)/df 2 < 0. The crucial assumption is that S c {f;T) vanishes at 
/min(r) with finite slope: 
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F[p(x)] 




Pa(x) P P (X) P(X) 



Fig. 7. Schematic illustration of the free energy landscape. For / m ; n < f < /ma* CO 
there are exponentially many local minima in an energy interval of width Af (two 
of them, p a (x) and pp{x), are shown explicitly). Above /max CO and below / m in(0 
there are at most an algebraic number of minima. 

^(/ raill ,T)<oc (78) 

(see FigurelHJ. At low temperatures (below T c ) the partition function Z(T, N) 
can be written as a sum over the individual local minima: 

Z(T, N) = exp[-/3JV#(T)] £* £ exp[-/3Nf a (T, JV)] (79) 



which becomes for N large 
Z 



(T,iV)« / d/exp{-JV[/3/-S c (/,r)]}& . (80) 



The main contribution to this integral comes from the minimum solution 
/* (T) of the free energy 

$(f,T)=f-TS c (f,T) , (81) 

i.e. 

$(T)=mm<?>(f,T) = r(T)-TS c (r(T),T) . (82) 
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The Boltzmann constant is put to one. Now, there are two possibilities. First, 
for temperatures below T c but high enough f*(T) will be within the interval 
[/min(T), / max (T)]. Then f*{T) follows from the solution of T)/8f = 0. 
Using Eq. ftHl this yields: 

^= d ~§fU,T) (83) 

(see Fig. EH). Second, f*(T) will decrease with decreasing T and will get stuck 
at f min (T) (see Fig. 133). Then it is S c {f*{T),T) = and 

= / mi „(T). (84) 

If we denote by So(T) the slope of dS c /df at / m i n (r) : 

So(T) = -^(/ min (T),T)<oc (85) 
this will happen at a temperature T s > with 

^r = so(T s ) . (86) 

Since S c vanishes at T s it is identical to the Kauzmann temperature Tr-, The 
reader should note that the existence of a nonzero Kauzmann temperature 
and accordingly the existence of a static glass transition at Tk relies on the 
finite value of sq(T) for all temperatures, including zero. 

Now, it is the major goal to calculate the free energy <P(T) and the config- 
urational entropy below Tk ■ This can be done by a trick. Instead of taking one 
system one chooses m replicas which are weakly coupled to each other with 
coupling constant e j3SHS3EI]- In the glass phase, i.e. below Tk this small 
coupling will force all systems into the same local minimum a. Therefore the 
corresponding partition function is given by: 

/mo»(T) 

Z m (T,N)^ J dfcxp{-N[mPf-S c (f,T)}} (87) 

/mi»(T) 

in analogy to Eq. JHqJ. The "saddle point" condition Eq. (183)) is: 

= = §(/, r> . (88) 

If we allow m to take any real positive value, it is obvious that Eci.f l88J) has 
a solution /„(T) > f m \ n (T) if m < 1 is small enough, even if /„ =1 (r) = 
/min(r). Then the free energy <P(m, T) is given by: 

*(m, T) = mm[m/ - T5 C (/, T)\ . (89) 
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Fig. 8. Schematic /-dependence of the configurational entropy S c (f,T). (a) Tk < 
T < T c where f*(T) is between / m i n (T) and f ma , x (T), (b) T < Tk where f*(T) 
coincides with f m m(T). Note that the slope of the tangent (which must equal 1/T) 
in (a) is smaller than in (b). The dashed line in part (a) shows the slope of the 
tangent at f m - m (T). 



Next, it is easy to prove that <£(m, T) and the free energy per particle of the 
replicated system: 



<j){m,T) = —<P(m,T) 



allows to calculate f(T) and S C (T): 

f(T) = 



d<P(m,T) 
dm 



(90) 



(91) 
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Sc(T) = ^ 9 -^ . (92) 

Hence, the knowledge of <P(m, T) allows to determine / and S c . Increasing to 
towards one there will be a critical value < m*(T) < 1 at which /^(T) = 
/min(r). A schematic representation of the m—T phase diagram is given in 
Fig. El There are two phases: a liquid phase above the solid line. For T > Tk 
it is a liquid where the replicas become uncorrelated for e — > whereas for 
T < Tk and to < m* (T) it is "molecular liquid" where the particles of the m 
replica form "molecules" even for e — > 0. Below the solid line there is a glass 
phase. Since S C (T) = in the glass phase, Eq. (192)) implies that <j)(m, T) is 
independent on to, i.e. it is: 

4>{m,T) = 4>{1,T) (93) 

for all to > to* (T) and T < Tk- On the other hand <j>(m,T) is continuous at 
the liquid-glass phase boundary to*(T): 

liquid (m* (T) , T) = glass (to* (T) , T) (94) 

combining Eq. Q93J1 and Eq. Q94J1 one arrives at the important result: 

${T) = 4>{l,T) = ^ icsai<i {rn*{T),T) . (95) 

Due to Eq. j^H one can calculate the free energy <P(T) of the physical system 
from the free energy of the replica system in its liquid phase, despite T < Tk- 
Since there are powerful techniques for the calculation of the liquid free energy 
|4()j the relationship Eq. II95H allows to calculate #(T) from <fi(m, T) and 
the latter also allows to determine the conflgurational entropy S c (f, T) from 
Eq. l|9ll and Eq. ll92l by eliminating to. 



4.2 Microscopic description 

In the last section we have shown that the free energy <P{T) and the conflg- 
urational entropy S c (f,T) can be obtained from <P(m,T), the free energy of 
a replica system in its liquid phase. Now we will describe how this can be 
performed from first principles. 

The potential energy of a N-particle system in a volume V may be given 
by pair potentials v(x): 

V(xi, xn) = ^ vjxj - Xj) . (96) 

i<j 

Let us consider to identical systems (replicas) which are weakly coupled by an 
attractive pair potential w(xf — Xj), which is considered to be short ranged. 
Then the total potential energy is given by: 
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Fig. 9. Phase diagram of the replica system in the m-T-space. Solid line separates 
the liquid from the glass phase (see also text). The dashed line within the liquid 
phase separates the "atomic" liquid for T > Tk — 0.5 from the "molecular" liquid 
for T < T K . The curved part of the solid line is m*(T). 



- E v (*? > ■ • • . *n) + £ E E w W 

o=l a<fc i<j 



»5) 



(97) 



where xf is the position of particle i in replica a and e > is an infinitesimal 
coupling constant. Note, that w breaks the replica permutational symmetry. 
It acts like a symmetry breaking magnetic field in case of a ferromagnet. For 
e > the thermodynamic limit N — > oo forces the replicas into the same local 
minima. Taking the limit e — > oo afterwards leaves the replicas in the same 
state for T < Tk and makes them uncorrected for T > Tk (cf. discussion 
in section 4.1). Therefore for T < Tk and e > it is (using an appropriate 
labelling of the particles) : 



Xi 



Xi 



Xi 



(98) 



Therefore we can introduce center of mass coordinates Xi of an "m-atomic 
molecule" and relative coordinates uf such that 



xt=X i + U t 



where 



^ m 



(99) 
(100) 



o=l 



and uf has to fulfil the N constraints: 
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E< = ° ( 101 ) 

a=l 

for all i. 

The classical partition function (configurational part) is given by: 

Z e m (T,N) = ^ J Y[d d xteM-^({xn)} (102) 

i,a 

where d is the spatial dimension. Making use of Eq. ((99]) this leads to 

Z s m (T, N) = ± ] Jl[ ^X, J J] d d u1 (103) 

i i.a 

d 

JJ (m^(^<)) exph/JV^tfX^K})] 

i— 1 a 

where the J-function accounts for the constraints Eq. H101J1 • Due to Eq. (l98l) , 
it is \uf\ <C 1. Therefore a harmonic approximation can be applied which 
yields: 

Z° m (T,N) = Z m (T,N) ^ / / I[d d u? ■ 

i i,a 

■ [] (m d «5(]T <) ) • exp { - /5[mF({X J) + 

i a 

+^EE - ^) 2 ] } (104) 



where (Mi W v) is the Hessian matrix of V and n? is the ^-component of uf. 
Note that we were allowed to put e = 0, because we already assumed that 
Eq. ll98l) holds. Using an integral representation of the (^-function the integral 
with respect to it? are Gaussian and can be performed which involves 



Nd 1/2 Nd 

[n **({*<})] [u x ^ x ^ 



k=l 



k=l 



-m/2 



exp 



m — 1 



IVlntM^^HJ,})) 



(105) 

where Afc are the eigenvalues of (Mi^j u ). The result can be brought into the 
following form: 



Z m {T,N) £* (c m (T)) Ar Z 1 (T/m,^)( exp 
where: 



to — 1 



t (T) =TO^(27rT)2 



|(m-l) 



Trln(Af VJI ,({XJ))J)(T/m) 
(106) 

(107) 
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Z\(T/m, N) is the partition function of the original system, but at a temper- 
ature T/m and ((•)) denotes averaging with respect to exp[— f3mV({Xi})]/Zi 
(T/m, N). In a next step one uses the approximation: 



exp 



exp 



-(TiHMi^dXi}))) . cms) 



The r.h.s. of Eq. 1110811 could be calculated if the density distribution of the 
eigenvalues would be known, which is not the case. Therefore, some more ap- 
proximations are necessary. We will not present these technical manipulations 
which can be found in Ref. jSQj but present the idea. Let us skip for a moment 
the logarithm in Eq. I|108jl . Then we have to calculate (Tr(M iM j V ({Xi}))) = 
£(M V , V ({X 4 })). Since = E^v/di^r = X { - X 3 ) this 

3 

involves the average of the second derivatives of the pair potential v(r). As- 
suming that Mi^in do not fluctuate much one obtains finally 

(Tr(M VJl/ ({X. i })))(^) « ^ J d d rg(r, ^) V 2 v{r) (109) 

i.e. (Tx{Mi^ y j V ({X i}))]) (T / m) can be related to the pair distribution function 
at T/m. In a similar way one can express (T^ln(Mi^j V ({Xi}))])(T/m) by 
g(r, T/m): 

(Trhi(M iw v({X i }))])(^) « Ng\g(r; ~)j . (110) 

\ / / V / L lib- 

Using these approximations we end up with the free energy <P{m,T) = 
— T lim iV"- 1 ln.Z m (T,j\n; 

0$(m, T) m - lnc m (T) + #(T/m) + i(m - l)0[ff(r; T/m)] . (Ill) 

This result demonstrates what has been discussed in section 4.1. The free 
energy <P(m, T) of the replica system is given by the free energy of the original 
system (m = 1) at T/m and by a functional Q of the pair distribution function 
of the original system, also at T/m. In the liquid phase <I>(T/m) can also be 
expressed by g(r 7 T/m) ^Hj- Therefore, up to an irrelevant term In c m (T) 
one has succeeded to express <P(m, T) by the pair distribution function at 
a temperature T/m. If T is below Tk, one has to choose m small enough 
in order to be in the liquid phase. As already mentioned above there are 
powerful methods 00] to calculate g(r, T/m) for the liquid phase. Finally, the 
such obtained ${m, T) allows to determine #(T) and S c (f, T) as described 
in section 4.1. In Figure E3 we show S c (f, T) obtained by Mezard and Parisi 
[49 for a three-dimensional system (d = 3) with density one and a soft sphere 
pair potential v(r) — | T" | 12 . The /—dependence has the qualitatively correct 
behavior, on which the discussion in section 4.1 has been founded. Using the 
condition Eq. l(8fijl for the static glass transition point with T s = Tk these 
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Fig. 10. Microscopic result for S c {f, T) for T = 0.05, 0.1, 0.15 (from left to right) 
for a system with soft sphere potential (see also text). 

authors obtained Tk — 0.194 or Tk — 1.51 for the dimensionless coupling 
parameter r = nT^ 1 / 4 . This value is in a reasonable range and is close to 
the values r^ m = 1.60 |£2] and 1.46 (£2] from a MD-simulation of binary 
soft sphere liquids. Since the authors of Ref. (231 claim that their results are 
compatible with MCT, it is not quite obvious whether the numerical values 
yield T c — nT c or not. 



5 SUMMARY 

In this article we have reviewed phenomenological and microscopic theories 
for the structural glass transition. The phenomenological approaches rely on 
several assumptions which are not proven to be correct. Although they are 
connected with appealing "physical pictures" their predictive power is limited. 
However, one class of them has obtained a microscopic justification by the 
replica theory for structural glasses |45l 1491 157)] . This theory, based on first 
principles, predicts a static glass transition at the Kauzmann temperature Tk 
where the configurational entropy S C (T) vanishes, as stated in the Adams- 
Gibbs-theory (1st chapter). Although one can avoid the use of replicas |5l] . 
the replica theory has a certain beauty because the several equivalent glass 
phases with same free energy can be described by copying the system m 
times and introducing a weak coupling between the copies (replicas). This 
coupling acts as a "symmetry" breaking field similar to a magnetic field for 
a ferromagnet. Analytically continuing m to positive real numbers allows to 
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relate the thermodynamical properties of the glass phase, i.e. for T < Tk, to 
those of the liquid phase, provided m is taken small enough with respect to 
one. 

Quite a different type of transition is obtained from mode coupling theory 
|27l 1281 I29j . MCT is a dynamical theory, in contrast to replica theory. It 
provides an equation of motion for, e.g. the spatial Fourier transform <£>(q, t) 
of the normalized density correlator of a simple liquid. Extension to molecular 
liquids is straightforward EDI EB E2 EBj ■ MCT predicts the existence of a 
dynamical glass transition at a temperature T c where the dynamics changes 
qualitatively. Above T c the correlator <P decays to zero, and converges to the 
nonergodicity parameter /(g) > 0, below T c . The nonergodicity parameters 
vary discontinuously at T c and can be interpreted as glass order parameters. 
Besides this, MCT makes several new predictions. Close to T c there exist 
two scaling laws, the critical and the von Schweidler law. The corresponding 
time scales t a (T) and r(T) exhibit power law divergence at T c . For times 
much larger than a typical microscopic time to and much smaller than the a- 
relaxation time scale t(T) the q- and i-dependence of $(q,t) factorize, which 
is a very strong statement. Furthermore, all the exponents of those power 
laws can be obtained from one parameter A, only, the exponent parameter. 
A is determined by the static correlator S(q) at T c . This proves the very 
microscopic nature of MCT. Although A and therefore the exponents are 
dependent on the physical system, i.e. on the interaction, they are universal 
for all correlators of one system which couple to the density fluctuations (21] ■ 

These two microscopic theories in some sense are complementary to each 
other and yield two different glass transition points. It is interesting that the 
existence of a static and a dynamical glass transition had also been found for 
mean field spin glasses with discontinuous order parameter. There it has even 
been speculated that spin glasses are quite similar to structural glasses [H]. 
For the mean field spin glasses both transitions are related to singularities 
and occur at Tk and T c with Tk < T c . Therefore both transitions are sharp. 
Concerning MCT (as described in the 3rd chapter), sometimes called idealized 
mode coupling theory, it has been shown [3U 1551 156] that the singularity is 
removed, due to ergodicity restoring processes. Nevertheless, if the time scale 
of these processes is much larger than t(T), one can observe in a certain 
time and temperature window the dynamical behavior as predicted from 
the idealized theory |27l 1281 I2fl| . Deviations start to emerge very close to 
T c . This fact has some importance for the static glass transition. For mean 
field spin glasses it was proven (S3 EHj that the system does not relax to 
equilibrium below T c . Therefore the static transition at Tk is masked. Due 
to the ergodicity restoring processes this is not anymore true for systems with 
finite range interactions. Despite of that, it is by far not clear that there is a 
minimum free energy f m \ n (T) at which the configurational entropy S c (f,T) 
vanishes linearly with a slope sq(T) — {dS c /df)(f m \ n (T),T) < oo. This is 
a very strong assumption which may not be fulfilled in general. For real 
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systems, i.e. in finite dimensions, with short range interactions the physical 
picture described in section 4.1 probably does not hold, i.e. there are not 
an exponential number of states with density p( a \x.) and infinite life time. 
However, there might exist systems with large enough frustration being close 
to the idealized situation, at least on a finite time scale. The same holds for 
MCT. For instance, the ergodicity restoring processes seem to be extremely 
weak for colloidal systems. Indeed it has been shown that the dynamics of 
colloids can be described over many decades in time by MCT [59| . 

Independent on whether the singularities of both microscopic theories ex- 
ist or not, the progress which has been made is significant, this is particularly 
true for MCT. The number of experiments and simulations |12l 1131 1141 115j 
(see also the contribution by U. Buchenau in this monograph) having found 
consistency with the MCT-predictions in an appropriately chosen time and 
temperature interval is enormous, despite the deviations very close to T c . One 
may hope that the few tests 60j of replica theory may be continued in order to 
check its validity in more detail. And an extension including time-dependence 
would be desirable as well. 

Both theories probably describe idealized situations only. As mentioned 
above MCT has been extended [HH EEj to include ergodicity restoring pro- 
cesses. But this extended version has not really been tested, because it is 
rather involved. Describing the ergodicity restoring processes by a single pa- 
rameter 8 there was a comparison with experimental data which is satisfac- 
tory ,61]. Whether replica theory can be extended in case that the singularity 
is spurious is not clear. Such an extension might also lead to a much more 
complicated mathematical structure. In that case it might be much better to 
restrict to the idealized version, which holds for both theories. 

Nevertheless there are still many challenging problems left. Let us mention 
some of them: 

(i) further investigations of similarities and dissimilarities between systems 
with and without quenched disorder 

(ii) glass transition in lattice gas models j£2] 

(iii) connection between the "potential energy landscape" properties j£[ E3 
IT%] with MCT and replica theory; e.g. it has been shown that the saddle 
index vanishes at the MCT-temperature T c 

(iv) investigation of models with trivial statics which may not exhibit a static 
glass transition but a MCT-like one |HH 

(v) investigations of models without or weak static correlations where 
replica theory and MCT (in its present form) do not yield a glass transi- 
tion, although found in simulations [66J 

(vi) nonequilibrium behavior (aging) 
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